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FOREWORD 

This  investigation  was  conducted  from  April  1,  1972  to  March 
30,  1975  by  the  Department  of  Electrical  Engineering,  University 
of  Minnesota,  under  Contract  DAAK02-72-C-0398  with  the  Advanced 
Research  Projects  Agency.  The  contracting  officers  representative 
is  Mr.  Edward  C.  Walker,  Night  Vision  Laboratory,  Fort  Belvior, 
Virginia.  The  final  technical  report  is  submitted  in  two  (2) 
parts. 


A.  van  der  Ziel 
Principal  Investigator 
J.  H.  Judy 
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Part  I 


A.  Theory  of  Pyroelectric  Detectors 
1.  The  Thermodynamical  Approach 

In  the  thermodynamical  approach  to  pyroelectricity  one  writes 
the  field  E  in  a  pyroelectric  as  an  odd  function  of  the  polariza¬ 
tion  P 

E  =  g(T  -  TC)P  +  BP3  +  CP5  +  ....  (1) 


where  Tc  is  the  Curie  temperature  and  8,  B  and  C  are  constants. 

We  assume  that  8  and  B  are  positive  and  for  the  sake  of  simplicity 
we  first  neglect  the  CP^term.  For  T  <  Tc  the  equation  E  =  0  has 
a  non-zero  solution  for  P;  in  other  words,  there  is  a  spontaneous 
polarization  Pg  for  T  <  Tc.  E  =  0  yields 

8  (T  -  T) 

8(T  -  TC)P  +  BP-3  =  0  ;  P  =  Ps  =  [ - - ]  '  (2) 


Hence  the  pyroelectric  coefficient  is 

_  3PS  _  1,6, 1/2  1 

'  3T  2  B  (Tc  -  T)l/2 

the  differential  susceptibility  is 


(3) 


e  -  1  = 


l/ef 


e03E/aP|p  8 (T  -  T_ )  +  3BP  2 


6e0(Tc  -  T) 


(1) 


and  the  figure  of  merit  of  the  pyroelectric  detector  is 

P  p  e01/2e 

e1//2  (e  -  1)1//2  (2B)1/2 

which  is  independent  of  the  temperature  T.  Experimentally  one 

1/2 

finds  some  decrease  in  p/(e  -  1)  for  decreasing  temperature  T; 

5 

most  of  it  comes  from  the  CP  term  in  (1). 


hi 


This  is  a  thermodynamical  theory/  and  hence  it  should  be 
correct.  Unfortunately  the  Taylor  expansion  (1)  does  not  con¬ 
verge  too  well  for  large  P,  and  as  a  consequence  (5)  (or  its 

5 

equivalent  if  the  CP  term  is  taken  into  account)  is  only  correct 
near  the  Curie  temperature.  Better  agreement  away  from  the  Curie 
temperature  might  be  obtained  by  using  still  higher  order  terms 
in  (1),  but  this  is  no  guarantee  that  the  Taylor  expansion  will 
converge  well  for  larger  P. 

To  carry  the  theory  further  it  is  necessary  to  express  the 

terms  6  and  B  in  molecular  parameters.  Only  then  can  one  hope 

1/2  -9  -2  -1 

to  understand  why  p/e  is  of  the  order  of  3  x  10  Coulomb  cm  °K 

We  shall  investigate  this  problem  in  the  next  two  sections. 

2.  Molecular  theories 

We  shall  now  discuss  two  theories  that  allow  the  expression 
of  6  and  B  in  terms  of  molecular  parameters.  To  that  end  one 
writes  for  the  potential  energy  V(x)  of  an  ion  in  the  crystal 

V (x)  =  ~  fQx2  +  i  gQx4  +  |  hQx6  +  ....  (6) 

where  f g,  g^  and  hQ  are  constants  and  x  is  the  deviation  from 
equilibrium  in  the  X-direction.  To  simplify  matters,  we  ignore 
the  sixth  order  term  in  (6).  We  then  have  two  possibilities: 

1)  f0  <  0,  gQ  >  0.  In  this  case  V(x)  has  a  minimum  when  dV/dx  =  0,  or 

fQx  +  gQx3  =  0  ;  x  =  xQ  =  +  -) 1//2  (7) 

so  that  the  ion  is  displaced  from  the  origin  and  thus  has  a  dipole 
moment.  If  the  ion  is  singly  charged,  the  dipole  moment  u  =  e|xQ|  = 
e(-fQ/gQ) 1//2  and  this  dipole  can  orient  itself  parallel  or  anti¬ 
parallel  to  the  X-axis.  We  shall  see  that  this  model  always  leads 


3 


to  ferroelec tricity.  It  is  called  the  dipole  model;  it  is  discuss¬ 
ed  in  the  next  section.  TGS  and  related  compounds  seem  to  be 
examples  of  it. 

(2)  fg  >  0,  >  0.  In  this  case  V(x)  has  its  minimum  value  at 

x  =  0,  so  that  the  origin  is  a  true  equilibrium  position  for  the 
ion.  When  the  interaction  with  the  other  ions  is  taken  into  account, 
it  may  be  shown  that  under  certain  conditions  a  spontaneous  dis¬ 
placement  of  the  ion  can  occur,  leading  again  to  ferroelectric 
behavior.  This  model  is  therefore  called  the  displacement- type 
model,  and  the  materials  satisfying  the  above  condition  are  called 
displacement-type  ferroelectrics .  They  are  discussed  in  Section  A. 4. 
Ba  Ti  0^,  Li  Ta  0^  and  SBN  seem  to  be  examples  of  it. 

3.  The  Dipole  Model 

We  assume  here  that  the  elementary  cells  in  the  material  can  be 
polarized  parallel  or  antiparallel  to  the  local  electric  field  E£ 
and  that  they  have  a  permanent  dipole  moment  p.  The  local  field 
acting  on  this  elementary  dipole  is 

Ep  =  E  +  M  (8) 

X,  Sg 

where  E  is  the  applied  field,  P  the  polarization  and  A  the  Lorentz 
factor. 

In  complete  analogy  with  the  magnetic  case,  it  can  be  shown 
that  the  net  polarization  P  for  N  elementary  dipoles  per  unit 
volume  is 

UEo 

P  -  Np  tanh{^)  (9) 

If  one  now  puts  E  =  0,  one  finds  that  there  is  a  temperature  Tc  so 
that  Eq.  (9)  has  a  non- zero-  solution  for  P  if  T  <  Tc#  In  other 
words  the  material  is  ferroelectric  for  T  <  T^ ,  and  Tc  is  the  Curie 
temperature  of  the  material. 
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We  now  invert  Eq.  (9)  and  substitute  for  E^.  This  yields 
after  a  Taylor  expansion  of  the  tanh  ^  term 


(E  +  ~)  =  P  +  ^  3 
^0  U  N  3p  r 


p J  + 

'  •  •  •  • 


kT, 


or  E  =  (*§-  -  |£)P  +  4  3 

M  N  b0  3y  NJ 


c  p3  + 

r  i  •  •  •  • 


(10) 


where  we  have  replaced  T  by  T^,  in  the  last  term.  Comparing  Eq.  (10) 
with  Eq.  (1)  we  see  that 


kT, 


3  = 


2 

M  N 


B  = 


o  4m3 

3m  N 


T  -  *  Q 
c  eo  k 


(10a) 


so  that 


(e  -  1) 


172  "  (2  e0  T 5 


Nkj 1/2 

C 


(11) 


We  now  see  that  the  density  of  dipoles  may  be  written 


N  = 


(12) 


where  AQ  is  Avogadro's  number,  p  is  the  density,  w  the  molecular 
weight,  and  n  is  the  number  of  molecules  per  unit  cell 
(this  takes  into  account  the  possibility  of  more  than  one  molecule 
per  unit  cell.)  This  yields 

- 2 — =  io.5  x  10-9( - £-= - )1/2  Coulomb  cm-2  °K~ 

(e  -  1)  n  JL  C 

n  100  1000 

-9  -2  -1 

This  usually  lies  between  (3-8)  x  10  Coulomb  cm  °K  whereas 

-9 

the  experimental  value  for  most  materials  lies  around  3  x  10 
Coulomb  cm  °K  .  For  TGS  near  the  Curie  temperature,  for  which 
the  model  should  especially  apply,  the  agreement  is  even  closer. 

In  view  of  the  rather  crude  model,  this  is  very  good  agreement. 
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This  theory  was  published  in  Ferroelec tries :  J.  H.  Judy, 

H.  D.  Park  and  A.  van  der  Ziel,  The  Molecular  Field  Model  Applied 
to  Pyroelectric  Detectors,  Ferroelectrics,  8,  685-687,  1974. 

4 .  The  Displacement-Type  Model 

We  start  with  the  potential  energy  of  the  ion 

V(x)  =  |  fQx2  +  i  gQx4  (14) 

The  equations  of  motion  of  the  ion  may  then  be  written,  if  the  ion 
is  assumed  to  be  singly  charged,  as 

"  3 

mx  +  fQx  +  gQx  =  e  E£  (15) 

where  m  is  the  mass  of  the  ion  and  E£  the  local  field. 

The  thermal,  incoherent,  anharmonic  oscillation  of  the  ions 

is  the  solution  of  the  homogeneous  equation 

"  3 

mx  +  f qX  +  gQx  =  0  (15a) 

with  the  solution 

x  =  Z  x„  cos  n  a)ftt  u)n  —  ——  (16) 

n  n  0  0  m 

We  must  now  solve  Eq.  (15)  under  the  presence  of  these  oscillations. 
To  that  end  we  write 

x  =  E  xfl  cOs  medgt  +  y^  (16) 

n 

and  average  over  a  complete  cycle.  We  then  obtain 

myl  +  [f0  +  390(Z  xn  cos  no)()t)2]y1  +  g^3  =  eE£ 
or 

H 

^1  +  ^1  +  ^1  =  eE£  (17) 
because  all  other  terms  disappear  in  the  averaging  process.  Here 
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3gnkT 

£  ■  fo +  V 


g  =  g. 


(17a) 


since,  according  to  the  equapartition  law 


(£  xn  cos  n  wQt) 
n 


2  „  kT 
■  ^ 


as  long  as  the  anharmonicity  is  not  too  large. 

We  shall  now  try  to  find  out  whether  there  is  a  case  for  which 
all  ions  are  displaced  over  the  same  distance  y^.  If  that  is  so, 
we  may  write 

vp 

(18) 


p  YeNy1 

E„  =  E  +  1—  =  E  +  - - 

1  eo  Eo 


Substituting  into  (17)  this  yields 

ah  +  (£  -  n  +  Vl*  - 


(19) 


If  a  temperature  Tc  can  be  found  such  that 

3gnkTr  e2wv 

t(Tc>  -  «o  +  ^ 


(19a) 


then  the  material  is  ferroelectric  for  T  <  Tc  and  paraelectric 
for  T  >  Tc#  The  necessary  and  sufficient  condition  for  this  is 


fn  <  Oi 

0  eo 


and  then 


TC  = 


(e  Ny/Eq  -  fQ)f0 
3g0k 


(19b) 


(19c) 


For  fp  >  e  Ny/Sp  the  material  is  always  paraelectric. 

If  condition  (19b)  is  satisfied,  Eq.  (19)  has  a  non-zero,  steady 
state  solution  for  E  =  0,  because 

(f  _  ejp^  +  gQyi3  =  o 


has  the  solution 
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*i  =  ^10  =  ( 


_  ,e  NY^£0  "  £0,l/2  _  ,3k(TC  '  T),l/2 

q„  '  ~  1  f  J 

yO  rO 


so  that  the  spontaneous  polarization  Pc  is 

D 

Ps  =  eNx10  =  eN(|^)1/2  (Tc  -  T)1/2 

and  the  pyroelectric  coefficient  p  is 

_  _  dPS  _  eN  ,3k, 1/2  1 

^  dT  2  f  ’  1 75 

0  (T  -  T)X/Z 

To  find  the  dielectric  constant  e,  we  apply  a  small  d.c. 
field  E,  so  that  y^  =  y^  +  Ay.  Neglecting  higher  order  terms 
Ay ,  Eq.  (19)  yields 


4y(f  +  S_SX  +  3g0  y102,=  eE 


or 


Ay  =  eE 


6g0K(Tc  -  T) 


and  AP  =  eNAy  =  e2N 


e2NT„ 


6g0k(Tc  -  T) 


E  = 


2(TC  -  T) (ezNy/e0  -  fQ) 


so  that  e  -  1  = 


AP  =  e_N  _ _ 

£0E  £0  2(TC  -  T) (e2NY/e0  -  fQ) 


and  uence  the  figure  of  merit  of  the  material  is 


(e  -  l)1'2  2Tc 


3e0N\l/2  ,e  NY/c0  "  f 0 » 1/2 
_  ) 


The  first  term  in  (22)  corresponds  to  the  expression  found  for 
dipole  model;  if  the  factor 


e  NvAq  -  fQ  1/2 
*  €  > 


(20) 

(20a) 

(20b) 

in 

(21) 

(21a) 

(21b) 

(22) 

the 

(22a) 


does  not  differ  too  much  from  unity  the  two  models  give  approximately 
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the  same  result.  It  again  makes  values  of  the  order  of  3.0  x  10 
-2  —1  1/2 

Coulomb  cm  °K  for  p/e  understandable. 

In  analogy  with  the  thermodynamical  model  one  might  expect 

that  higher  order  terms  in  V{x)  would  give  rise  to  a  significant 

1/2 

temperature  dependence  of  p/e  .  However  this  has  not  been  worked 
out  in  detail. 

This  paper  is  being  published  in  one  of  the  forthcoming 
issues  of  Ferroelectrics : 

1/2 

H.  D.  Park,  A.  van  der  Ziel  and  J.  H.  Judy,  p/ (e  -  1) 
for  Displacement  Type  Ferroelectrics,  Ferroelectrics,  in 
the  press. 

The  work  contained  in  Section  A  has  appeared  in  greater 
detail  in  Mr.  H.  D.  Park's  Ph.D.  thesis. 
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B .  Other  Theoretical  and  Experime ntal  Work  on  Ferroelectrics 
and  Dielectrics 

1 .  Feedback  Theory  of  Ferroelectricity 

We  establish  here  an  analogy  between  the  theory  of  ferro¬ 
electricity  and  feedback  theory. 

Let  there  be  N  polarizable  atoms  per  unit  volume  in  a  dielectric, 
each  having  a  polarizability  a,  then  the  polarization  P  is 


P  =  NaE0  =  Na  (E  +  ^-)  ; 


(23) 


where  is  the  local  field,  E  the  applied  field  and  y  the  Lorentz 
factor.  Solving  for  P  yields  for  the  susceptibility 

r,  Na/e. 


e  -  1  = 


eQE  1  -  yNa/ eQ 


(24) 


corresponding  to  the  gain  of  a  feedback  amplifier  with  a  gain 
Nct/£Q  for  zero  feedback  and  feedback  factor  y. 

If  yNa/e^  <  1  at  all  temperatures,  the  material  is  a  normal 
dielectric;  if  yNa /Eg  >  1  for  some  temperature  range,  then 
the  material  has  spontaneous  polarization,  or  is  ferroelectric, 
in  that  temperature  range.  These  criteria  correspond  to  the 
Nyquist  condition  in  feedback  amplifiers. 

In  the  ferroelectric  regime  there  is  polarization  at  zero 
applied  field.  According  to  (23)  this  is  only  possible  if  u  is 
a  function  of  P  as  well  as  of  the  absolute  temperature  T,  such  that 


yNa (T,P)  E  x  (2 

e0 

in  that  temperature  range.  This  means  that  a(T,P)  decreases  with 

increasing  |P|  and  that  |p|  adjust  itself  such  that  (25)  is  sat- 

2 

isfied;  it  also  means  that  a(T,P)  must  be  a  function  of  P  . 


Equation 


(25)  is  equivalent  to  the  Nyquist  condition  for  a  stable  amplitude 
in  oscillators. 


For  most  ferroelectric  materials  there  is  spontaneous  polariza¬ 
tion  Pg  below  the  Curie  temperature  Tc  and  Pg  =  0  at  T  =  Tc.  The 
Curie  temperature  is  then  given  by 


YNaQ 


e 


(T  ) 

— _  =  i 
0 


(26) 


when  ctg  (T^)  is  the  value  of  a(T,P)  for  P  =  0.  Corresponding  to 
the  limit  of  stability  in  feedback  amplifiers. 

In  analogy  with  the  feedback  amplifier  we  may  write 

(T) 

a(T,P)  =  - 2__ - t -  (27) 

1  +  BP  +  CP  +  .... 


Substituting  into  (23)  and  solving  for  the  field  E  yields 

B  p3  +  C 


1  -  yNan(T)/en 
E  =  - ..  - -  P  + 


NaQ (T) 


NaQ(T) 


NaQ(T) 


P5  + 


(27a) 


which  corresponds  to  the  thermodynamical  theory  of  ferroelectricity . 
That  theory  thus  corresponds  to  what  would  be  expected  from  feedback 
considerations . 

Since  the  thermodynamical  theory  has  now  been  established, 

one  can  calculate  e  -  1  and  the  pyroelectric  coefficient  p,  near 

1/2 

the  Curie  temperature  and  one  can  evaluate  p/(e  -  1)  .  It  is  also 

easy  to  show  that  ctg(T)  must  decrease  with  increasing  T  for  T  near 

V 

This  work  is  scheduled  to  be  published  in  one  of  the  forthcoming 
issues  of  Ferroelectrics : 

A.  van  der  Ziel,  J.  H.  Judy  and  H.  D.  Park,  Feedback  Theory 
of  Ferroelectricity,  Ferroelectrics  9,  1975,  in  the  press. 


h 
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2 .  Solar  Power  Generation  with  the  Pyroelectric  Effect 
This  work  was  performed  while  at  the  University  of  Florida 
and  only  the  publication  costs  were  charged  to  the  contract. 

As  is  well  known,  when  an  incoming  power  density  P  is  chopped 
at  the  rate  f/sec  and  exp(jajx)  is  the  power  density  of  frequency  w 
then  the  current  amplitude  is 

nApP, 

Id0  =  cd (1  +  i/jwr „)  (28) 

ri 


where  to  =  2irf,  n  is  the  emissivity  of  the  front  face,  p  the  dielectric 
coefficient,  c  the  specific  heat  capacity  per  unit  volume,  A  the 
device  area,  d  the  device  thickness  and  t„  the  heat  time  constant 

n 


of  the  system 


TH  = 


c  d 


(n  +  n'MaT' 


(28a) 


where  n*  is  the  emissivity  of  the  back  face,  T  the  absolute  temper¬ 
ature  of  the  device  and  a  the  Stefan-Boltzmann  constant. 

The  device  has  a  capacitance  C  =  eeQA/d.  If  a  load  resistance 
R  is  applied,  maximum  power  is  transferred  if  wCR  =  1  and  the  trans¬ 
ferred  power  is 
P 


do 


n2 P2 (p2a)2 


max 


4toC 


4c2d2 [1  + 


.  WEE  A 

1  ]  — 


2  2 
U)  th 


which  has  an  optimum  value  P  .  if  wt,,  =  1,  so  that 

opt  H 


opt 

A 


n2P2P, 


8c  (n  +  n * )  4aT  ee, 


P  ./A 

The  efficiency  n  =  — *-«■ — 

Jr  r 


(29) 


(30) 


turns  out  to  have  a  very  small  value.  This  is  not  so  surprising 
for  the  device  is  essentially  high  impedance  device;  it  gives 


a  large  voltage  but  only  a  small  power.  Nevertheless  it  was 
thought  worth  while  to  put  it  on  record;  to  our  surprise  it  gained 
a  lot  of  interest. 

This  work  was  published  in  the  Journal  of  Applied  Physics: 

A.  van  der  Ziel,  Solar  Power  Generation  by  the  Pyroelectric 
Effect,  J.  Appl.  Phys.,  4_5,  4128,  Sept.  1974. 


3.  Equivalent  Circuit  and  Equipartition  Theorem  in  Dielectrics 
At  high  frequencies  the  dielectric  constant  e  becomes  complex 
and  its  absolute  value  is  less  than  the  low-frequency  value  e  , 

It  was  shown  that  for  a  dielectric  without  low-frequency  losses 
the  susceptibility  may  be  written 

e  -  1 


e  -  1  = 


1  +  jlOT 


(31) 


where  x  is  the  time  constant  of  the  dielectric,  so  that  the  ad¬ 
mittance  of  the  capacitor  is 


Y  = 


^eeQ  d  " 


A 

3O)E0  J  + 


ju>e0(es  -  1 ) A 
d(l  +  ju)x) 


(32) 


This  leads  to  the  equivalent  circuit  of  Fig.  la,  with 

r,  _  Td 


C0  = 


e0A 


e0(Es  "  1)A 


ps 


ps 


e0(es  "  1)A 


(32a) 


The  resistance  R  shows,  of  course,  thermal  noise, 
ps 

Calculating  the  mean  square  output  voltage  yields 


£g 

+  c 


0 


(33) 


indicating  that  the  equipartition  law  does  not  hold  in  this  case. 

This  is  not,  however,  what  is  measured  when  the  capacitor  is 
a  pyroelectric  capacitor  that  is  part  of  an  array  of  pyroelectric 


detectors.  Let  this  array  be  scanned  at  the  rate  fg,  then  the 

lower  cut-off  frequency  of  the  system  f^  is  of  the  order  of  fg 

and  the  upper  cut-off  frequency  f ^  is  determined  by  the  scanning 

time  of  a  single  element  of  the  array.  We  must  now  evaluate  the 

2 

contribution  dAVg  to  the  output  voltage  between  the  frequencies 
f  and  f  +  df  and  integrate  over  the  passband  of  the  system 


AVq2  =  /  2  d(AV02) 

tl 


(33a) 


Usually  <<  1  and  then  AV^  is  much  smaller  than  kT/Cs. 

This  comes  about  because  in  the  previous  case  the  main  contribution 
2 

to  AVq  came  from  the  very  high  frequencies  to  which  the  system 
does  not  respond. 

For  a  dielectric  with  low-frequency  losses  it  was  shown  that 
the  susceptibility  may  be  written 


=  £°°  ~  1  £s  ~  > 

L  1  +  jtiJT  1  +  jo>T, 


( >>  T) 


where  e  -  1  is  the  low  frequency  susceptibility  (cox..  <<  1)  and 

b  -L 

e  -  1  is  the  high-frequency  susceptibility  (1/x.  <<  go  <<  1/x)  . 
This  expression  is  reasonably  accurate  if  the  low  frequency  losses 
with  the  relatively  long  time  constant  x^  are  not  too  large 
(e  -  e  <<  e  -  1). 

g  oo  oo 

The  device  admittance  is  in  this  case 


j«e0A  _  -  1)eqA  _  jw(eg  -  ^JcQA 

Y  -  d  +  d(l  +  j (a) x )  +  d(l  +  jcox^) 


This  can  be  represented  by  the  equivalent  circuit  of  Fig.  lb  where 


= 


£0(£~  ~  1)A  .  _  _  Td 

d  '  P°°  £0(£c=  "  1)A 


en<es  "  eco)A 


(35a) 
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The  resistors  and  R  show,  of  course,  thermal  noise. 

2 

If  one  calculates  the  contribution  to  AVQ  due  to  the  thermal 
noise  of  R^,  and  one  neglects  R  ,  one  obtains 


2  kT  £s  £° 


AV0  = 

s 


(36) 


This  is  much  less  than  the  equipartition  law  would  indicate,  and 
the  contribution  of  the  dielectric  losses  to  the  output  voltage  in 
a  pyroelectric  array  would  be  correspondingly  small. 

We  can  also  put  the  latter  as  follows.  Normally  one  char¬ 
acterizes  the  device  by  its  capacitance  C  and  the  losses  by  the 
loss  tangent  tan6.  Over  the  limited  range  of  frequencies  usually 
used  in  arrays  C  and  tan6  are  practically  constant.  Therefore  the 
open  circuit  device  voltage  has  a  spectrum 


Sv(f)  =  4kT 


tan6 

toC 


(37) 


so  that  the  mean  square  output  voltage  is 


AV 


ef  f 


2  =  /  2Sv(f)df  =  |  ~  Hn(^-)tan6 

fl  1 


(37a) 


This  can  be  much  less  than  kT/C,  the  value  given  by  the  equiparti¬ 
tion  law. 

This  work  was  published  in  the  Journal  of  Applied  Physics: 

A.  van  der  Ziel,  Equivalent  Circuit  and  Equipartition  Theorem 
in  Ideal  Dielectric  and  Ferroelectric  Capacitors,  J.  A.  P.  45 , 
1400-1401,  March  1973. 

A.  van  der  Ziel,  Equivalent  Circuit  and  Equipartition  Theorem 
in  Lossy  Dielectric  and  Ferroelectric  Capacitors,  J.  A.  P.  45, 
1402-1403,  March  1973. 
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4 .  Partition  Noise  as  the  Cause  of  Thermal  Noise  in  a  Two- 
Level  Ferroelectric 

In  a  two-level  ferroelectric  the  dipole  moment  can  orient 
itself  parallel  or  antiparallel  to  the  local  electric  field  E£  = 
E  +  yP/Eq,  E  being  the  applied  field,  P  the  polarization  and  y 
the  Lorentz  factor.  Let  N  be  the  density  of  dipoles  and  let 
N^  be  parallel  to  the  field  and  N2  be  antiparallel  to  the  field 
so  that  N^  +  N2  =  N.  The  polarization 


P  =  W  (n2  -  n2) 


AP  =  2yAN, 


(38) 


since  AN2  =  -  AN^.  Let  N^g  and  N2g  be  the  equilibrium  values  of 
N-^  and  N2. 

Now  according  to  the  theory  of  polarization  noise  the  polariza 
tion  noise  spectrum  is  given  by 

.2 


SAP(f>  = 


4kTeQ (eg  -  1) 
Vt 


.  .  2  2 
1  +  0)  T 


(39) 


(Nyquist  equivalent)  where  is  the  low  frequency  dielectric 
constant,  t  the  macroscopic  time  constant  of  the  dielectric  and 
V  the  volume  of  the  sample.  Working  backwards,  it  is  shown  that 


the  fluctuation  AN^,  without  the  feedback  caused  by  the  dependence 
of  the  local  field  upon  the  polarization,  is  given  by 

- =■  N.  »N, 


AN, 


’10‘  20 


N10  +  N20  V 


(40) 


which  is  just  the  partition  noise  expected  for  the  two-level  system 
The  conclusion  is  therefore  that  the  thermal  noise  in  this  case 
is  caused  by  the  partition  noise  of  the  two  level  system. 

This  work  was  published  in  Ferroelectrics : 

A.  van  der  Ziel,  H.  D.  Park  and  J.  H.  Judy,  Partition  Noise 
as  the  Cause  of  Thermal  Noise  in  a  Two-Level  Ferroelectric, 
Ferroelectrics,  8,  689-690,  1974. 
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5.  Polarization  Fluctuations  in  Capacitive  Bolometers 
A  capacitxve  bolometer  uses  a  ferroelectric  operating  slightly 
above  the  Curie  temperature.  Modulated  radiation  incident  upon 
the  bolometer  gives  rise  to  a  fluctuating  temperature  and  hence  a 
fluctuating  capacitance,  and  this  fluctuating  capacitance  is 
detected  by  applying  a  d.c.  voltage  V  to  the  capacitor,  corres¬ 
ponding  to  an  applied  field  E  =  V/d. 

At  zero  field  the  noise  resistance  R^  of  the  open-circuit  cap¬ 
acitance  C  was  equal  to  the  series  resistance  Rs  =  tan6/ (wC)  where 
co  is  the  angular  frequency  and  tan6  the  loss  tangent  of  the  material, 

but  at  high  fields  R  >  R^  and  R  saturated  at  high  fields.  There 
3  n  s  n 

is  presently  no  satisfactory  explanation  for  this  phenomenon. 

This  work  was  published  in  the  Journal  of  Applied  Physics: 

V.  P.  Singh  and  A.  van  der  Ziel,  Polarization  Fluctuations 
in  Capacitive  Bolometers,  J.  A.  P.  £5,  1452-1453,  March  1974. 


6.  Limiting  Flicker  Noise  in  MOSFETs 

It  is  very  tempting  to  use  MOSFETs  in  pyroelectric  arrays, 
since  they  are  more  easily  integrated  into  the  system.  Single 
pyroelectric  detectors  usually  use  JFETs  since  they  give  lower 
noise.  It  is  therefore  important  to  point  out  that  MOSFETs  are 
inherently  noisier  than  JFETs. 

The  reason  is  a  very  simple  one.  The  gate-oxide-semiconductor 
capacitance  C  has  a  loss  tangent  tan5  and  hence  a  series  resistance 
R s  =  tan6/(u>C)  that  shows  thermal  noise.  The  limiting  noise  re¬ 
sistance  R„  of  a  MOSFET  is  therefore 


Rn  =  *s  = 


tan  6 

GdC 


(41) 


-4  6 

If  we  take  tan6  =  10  and  C  =  5pF,  then  R^  =  3  x  10  /f  ohm. 


corresponding  to  3Mft  at  1  Hz.  In  contrast,  good  JFETs  may  have 
Rn  =  104fi  at  1  Hz. 

This  work  is  being  published  in  Solid  State  Electronics: 

A.  van  der  Ziel,  Limiting  Flicker  Noise  in  MOSFETs ,  Solid 
State  Electronics,  in  the  press. 
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C.  The  Capacitive  Bolometer 

1 .  The  D.C.  Capacitive  Bolometer 

For  a  ferroelectric  capacitor  operating  slightly  above  the 
Curie  temperature  T^,,  the  material  is  paraelectric  and  the  cap¬ 
acitance  is 


const 
T  -  Tc 


(42) 


Modulated  radiation  of  frequency  ui  produces  a  variation  in  cap¬ 
acitance,  this  variation  is  detected  by  applying  a  d.c.  voltage 
Vq  in  series  with  the  capacitor  and  connecting  the  circuit  to  a 
large  load  resistance  R  such  that  to  CqR  >>  1,  where  is  the  equil¬ 
ibrium  value  of  C. 

The  heat  response  of  the  detector  due  to  the  incident  radia¬ 
tion  exp(ju)t)  is  given  by 


'H 


dAT 

dt 


+  (9h  +  )AT  =  nP,  exp(jwt) 


H 


(43) 


where  the  heat  capacitance  CH  =  cAd,  where  c  is  the  specific  heat 
per  cm^,  A  the  device  area  and  d  the  device  thickness,  n  is  the 
emissivity  of  the  front  face,  g„  and  g  '  are  the  heat  loss  con- 

n  n 

ductances  of  the  front  and  back  face  and  AT  is  the  temperature 
change  due  to  the  incident  radiation. 

Putting  AT  =  ATq  exp(jcot)  and  solving  for  ATQ  yields 

nP, 


AT0  "  jtt>CHll  +  l/(  jutH)  ] 


(44) 


where  =  q  / (gH  +  g  ' )  is  the  thermal  time  constant  of  the  system. 

Since  =  1  sec  and  to  =  10  -  100  per  sec.,  to^T^  >>  1,  so  that 


AT, 


"1 

JtoC 


H 


ju)cAd 


(44a) 


The  capacitance  may  therefore  be  written 
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c  =  co  +  aT^  ATo  exP<:iut>  (45) 

If  the  capacitance  is  reasonably  linear  with  respect  to  the  applied 


voltage,  the  charge  Q  of  the  capacitor  is  CVg  and  the  short-circuit 
current  is 

gc 

1  =  3t  =  V0  51  =  jajV0  TT  AT0  exp(jwt)  (46) 

Hence  the  output  voltage  Vq  exp(jwt)  has  an  amplitude 


_  jalV0OC0/3T)AT0  3^  ^1E0 

0  JmCT  '  v0  C0  3T  *0  jiucA  (T  -  Tc) 


where  Eg  =  Vg/d ,  so  that  the  output  voltage  is  proportional  to 
the  d.c.  field  Eg .  For  large  field  saturation  effects  will  set 
in  that  will  be  discussed  in  a  moment. 

Assuming  the  noise  of  the  device  to  be  thermal  noise  of  the 
dielectric  losses 

c  4kT  tan6  4kT  tan6  d  ,, 

Vm  u)C0  "  u>eeQA  1 

Defining  the  noise  equivalent  power  Pe<^  by 


(Vrms  =  [SV(f>> 


yields 


P  .  C(^~  (4%tan6)l/2(ujAa)l/2  (49) 

eq  nE0  ccQ 

-1/2  3/2 

Since  e  varies  as  (T  -  T  )  ,  P  varies  as  (T  -  T  )  .  The 

c  eq  c 

lowest  Pecj  is  thus  obtained  by  keeping  the  temperature  as  close 
to  the  Curie  temperature  as  possible. 

When  saturation  effects  set  in,  one  must  use  the  charge  Q  =  AP(VQ,T) 
directly,  where  P(Vg,T)  is  the  induced  polarization.  The  current 
is  therefore 


I(t)  =  I.n  exp(jwt)  =  jioA  4=  AT.  exp (jut) 


(50) 
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and  the  output  voltage  has  an  amplitude 


v  -  dO 

0  JldC 


=  A_  3P  ,  1_  3Q  AT 
0  CQ  3T  0  CQ  ST  0 


where 


(51) 


C 


0 


(52) 


Q  = 


We  now  introduce  the  d.c 
CdV0*  Hence 


v 

0 


capacitance 


3Cd 

ST 


AT 


0 


Q/VQ,  so  that 


(53) 


Eq.  (47)  must  now  be  replaced  by  (53),  i.e.  3Cq/3T  must  be  replaced 

by  SC^/ST.  Saturation  means  that  the  factor 

i  3C, 

V  i—  — — 

0  CQ  3T 


saturates. 

At  large  fields  the  noise  is  somewhat  larger  than  the  thermal 
noise  of  the  dielectric  losses.  pegf  as  given  by  (49),  is  thus 
somewhat  too  optimistic. 

Experimental  results  are  discussed  in  Section  C3. 


2 .  The  A.C.  Capacitive  Bolometer 

In  the  a.c.  biased  bolometer  an  a.c.  signal  Vq  cos  w^t  is 
applied  to  the  bolometer,  where  >>  the  frequency  u>  of  the  in¬ 
cident  radiation,  and  series  tunes  the  circuit  by  an  inductance  L 

connected  between  the  a.c.  signal  and  the  bolometer.  If  Cq  is 

2 

the  quiescent  capacitance  of  the  bolometer  this  means  that  u>q  LCq  =  1. 
The  output  signal  is  taken  from  the  capacitor. 

We  can  now  calculate  the  response  and  find  that  it  is  the 
same  as  for  the  d.c.  case,  except  that  the  incident  radiation 
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produces  two  sidebands  of  frequencies  uig  -  w  and  Wg  +  w,  that  can 
be  detected  by  a  linear  detector.  The  noise  emf  per  unit  band¬ 
width  is  series  with  the  capacitance  now  has  a  spectrum 


Vf>  "  <4kT 


(54) 


Consequently  if  one  calculates  P  ,  one  finds 


(P  )  =  (P  ),  (——) 1//2 

eq  a.c.  eq  d.c.  ajg 


(55) 


For  u/co  =  1000  the  a.c.  value  of  P  is  thus  30  times  smaller 

eq 

than  the  d.c.  value,  indicating  the  great  advantage  of  the  a.c. 
bolometer. 

This  scheme  has  not  been  tried  in  practice.  There  are  some 
difficulties  to  be  overcome;  for  example  the  capacitance  C  not 
only  varies  in  the  rhythm  of  the  signal  of  frequency  w  but  also  in 
the  rhythm  of  the  signal  of  frequency  It  would  be  worth  while, 

however,  to  try  it  out  in  practice  and  see  what  can  be  achieved. 


3 .  Experimental  Results  for  the  D.C.  Bolometer 

In  his  Ph.D.  thesis.  Dr.  Singh  has  given  extensive  measurements 
of  the  capacitive  bolometer  effect  in  TGS.  He  finds  that  the 
sensitivity  in  V/W  is  linear  to  the  d.c.  field  EQ  at  low  fields  but 
saturates  for  fields  higher  than  lOkV.  Closer  to  the  Curie  temper¬ 
ature  saturation  effects  set  in  earlier.  He  also  finds  that  the 
inverse  of  the  sensitivity  varies  linearly  with  T  -  T^,  as  expected 
from  Eq.  (47)  .  He  has  also  measured  the  effect  of  the  d.c.  field 
on  the  small-signal  capacitance.  He  compared  the  series  resistance 

R_  and  the  noise  resistance  R  of  the  device  and  found  >  R_  at 
s  n  n  s 


high  d.c.  fields 
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Finally,  he  has  compared  the  operation  of  sample  TGS  42,  which 

had  a  Curie  temperature  of  49.8°C,  as  a  capacitive  bolometer  at 

52.7°C  and  as  a  pyroelectric  detector  at  room  temperature.  He 

found  in  the  first  case  a  maximum  noise  equivalent  however  P  = 

-9  1/2 

1.26  x  10  W/Hz  '  at  a  field  of  9kV,  whereas  the  same  device, 

when  used  as  a  pyroelectric  detector  at  room  temperature  had  P  = 

-9  1/2 

1.71  x  10  W/Hz  .  Since  operation  closer  to  the  Curie  temperature 
would  have  been  feasible,  and  this  would  have  resulted  in  a  lower 
Peq,  it  may  be  concluded  that  the  capacitive  bolometer  is  intrin¬ 
sically  able  to  compete  with  the  pyroelectric  detector. 

-9  1  /2 

Dr.  Singh  also  compared  the  value  of  P  =  1.26  x  10  W/Hz  ' 

with  the  theory  as  presented  by  Eq.  (49) .  He  obtained  a  theoretical 

value  for  P  that  was  about  a  factor  8  smaller.  This  factor  comes 
eq 

from  two  reasons 

a)  Saturation  effects  reduce  the  sensitivity  at  high  fields 
from  the  value  calculated  by  Eq.  (47). 

b)  At  high  fields  the  noise  resistance  Rn  is  considerably 

larger  than  the  series  resistance  R&  of  the  device,  and  this  means 

that  Eq.  (49)  must  be  multiplied  by  the  factor  (R  /R  )1,/2. 

n  s 

Both  factors  together  can  probably  explain  the  data. 

The  details  of  this  work  are  presented  in  the  following  thesis: 

Vijay  Pal  Singh,  Noise  in  Ferroelectric  Detectors  for  Infrared 
Radiation,  E.  E.  Department,  University  of  Minnesota,  December  1973. 
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D.  Other  Work  Performed  Under  the  Contract 

1 .  Diffusion  Theory  of  the  Response  of  Pyroelectric  Detectors 
In  a  series  of  two  papers  the  response  of  pyroelectric  detectors 
was  discussed  as  an  application  of  the  diffusion  theory  for  heat 
flow.  The  response  was  expressed  in  terms  of  hyperbolic  functions, 
which  greatly  simplifies  the  notation  and  makes  it  possible  to 
understand  the  frequency  response  almost  by  inspection. 

In  the  first  paper,  co-authored  with  S.  T.  Liu,  the  theory 
was  applied  to  a  single  free-bearing  detector.  At  very  low  and 
very  high  frequencies  the  results  agree  with  the  elementary  lumped 
circuit  approach  to  the  heat  problem;  at  intermediate  frequencies 
the  new  result  improves  on  the  lumped  circuit  model. 

In  the  second  paper  the  theory  was  applied  to  thin  pyro¬ 
electric  films  on  a  substrate  and  the  response  was  calculated  as 
a  function  of  the  substrate  and  film  thickness.  It  was  found  that 
for  low-frequency  applications,  the  pyroelectric  layer  should  not 
be  made  too  thin,  preferably  several  micrometers  thick,  that  the 

substrate  should  be  as  thin  as  possible  and  that  the  product  of 

3 

its  specific  heat  per  gram,  its  density  in  grams/cm  and  its  heat 
conductivity  should  be  as  small  as  possible.  For  thick  substrates 
there  is  an  intermediate  frequency  region  where  the  value  of  D*  is 
practically  independent  of  frequency.  This  intermediate  frequency 
region  extends  to  higher  frequencies  if  the  pyroelectric  layer  is 
made  thinner,  but  D*  decreases  with  decreasing  film  thickness 
in  that  region.  The  second  study  was  made  in  support  of  our  sputter¬ 
ed  thin  film  program.  These  papers  were  published  in  the  Journal 
of  Applied  Physics; 

1.  A.  van  der  Ziel  and  S.  T.  Liu,  Diffusion  Theory  of  the 

Response  of  Pyroelectric  Detectors,  J.  A.  P.,  £3#  4260,  1972. 
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2.  A.  van  der  Ziel,  The  Pyroelectric  Response  and  D*  of  Thin 
Pyroelectric  Films  on  a  Substrate,  J.  A.  P.  A5_,  546-549, 
1973. 

2.  Thin  Film  Program  (Dr.  Wehner) 

Films  (0.1  -  8 u  meter  thick)  of  lead-zirconium-ti tanate  oxide 
were  prepared  by  RF  sputtering  in  pure  Ar  and  a  10%  0^  in  Ar  mixture. 
Substrates  of  glass,  glazed  alumina,  mica  and  fused  quartz  were 
heated  to  500°C  during  deposition.  The  film  deposited  in  pure  Ar 
became  dark  grey  whereas  those  prepared  in  the  10%  mixture  were 
yellow  but  quite  transparent.  Capacitance  measurements  versus 
temperature  indicated  a  Curie  temperature  of  about  160°C.  Auger 
spectroscopy  studies  indicated  that  the  thin  film  had  less  Pb 
than  the  bulk  (about  30%). 

PLZT  (8%  La,  65%  Zr,  35%  Ti)  were  sputtered  as  thin  layers 
(0.1  -  0.5  y  meter  thick)  on  fused  quartz  and  mica  substrates  in 
a  pure  Argon  RF  triode  system.  The  substrate  temperature  was  varied 
from  20°C  to  500°C.  Auger  spectroscopy  studies  indicated  that  the 
films  were  deficient  in  La  and  Zr  to  the  extent  of  about  45%  less 
than  the  target  composition.  The  films  appeared  to  be  ferroelectric, 
as  capacitance  measurements  indicated. 

The  studies  indicated  in  general  that  while  it  is  possible 
to  deposit  layers  that  are  ferroelectric,  it  is  not  so  easy  to 
deposit  them  with  the  same  composition  of  the  bulk.  While  such  a 
deposition  would  be  worth  while,  it  will  probably  require  a  pro¬ 
longed  effort  to  make  acceptable  layers. 

3.  Materials  Program  (Dr.  Wertz) 

The  materials  program  produced  a  variety  of  crystals  useable 
as  pyroelectric  detectors.  The  most  useful  crystals  were  grown 
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crystals  of  TGS;  crystals  of  rather  good  quality  were  obtained  in 
this  manner.  Later  TGS  crystals  were  bought  from  Isomet. 

We  also  made  crystals  from  pressing  powders  at  high  pressure 
and  temperature  by  the  Honeywell  Corporation.  Ba  Ti  O^and  (Ba,Sr)Ti  0^ 
crystals  were  a  good  example  of  such  a  procedure.  Other  crystals 
were  Ba (Sn,Ti) 0^ /  with  a  range  of  compositions  of  the  vicinity  of 

Ba  (Tias'^o.lS*  °3* 

We  also  made  crystals  of  Tl^  As  Se^»  which  has  a  relatively 
small  gap  width,  and  is  a  semiconductor  rather  than  an  insulator, 
and  the  related  compound  In^  As  Se^,  which  was  expected  to  have  a 
larger  gap  width  and  would  be  better  at  room  temperature.  None 
of  the  crystals  were  large  enough  to  make  good  detectors. 

We  also  prepared  Pb,.  (Ge, Si)  ^0^  which  has  a  Curie  temperature 
of  84°C.  Only  small  crystals  were  made  by  the  Crochzalski  technique. 

Later  PLZT  and  SBN  samples  became  available  from  Honeywell 
Corporation  and  other  sources  and  many  devices  were  made  of  these 
materials. 

4 .  Contact  Problems 

At  first  we  had  some  serious  problems  with  contacts.  Especially 
with  TGS  it  was  very  difficult  to  make  good  devices  with  silver 
or  platinum  paste.  Apparently  the  solvent  attacked  the  crystal, 
deteriorated  its  performance  and  made  the  device  unstable. 

Later  very  satisfactory  results  were  obtained  by  vacuum  deposition 
(pressure  lower  than  10"_6Torr).  First  a  thin  layer  of  chronium  was 
evaporated  onto  one  side  of  the  device,  followed  by  a  thicker 
layer  of  gold.  Then  the  device  was  turned  over  and  the  process 
repeated  on  the  other  side.  Thin  gold  wires  were  (Sigmund  Cohn 
wire,  1-P)  attached  to  the  two  electrodes  with  the  help  of  silver 
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epoxy  (Dupont,  Silver  composition  4922)  and  the  device  was  mounted. 
This  produced  good,  stable  devices. 

Samples  that  could  withstand  a  much  higher  temperature,  like 
SBN  and  PLZT,  were  first  contacted  with  the  help  of  burned  silver 
or  burned  gold.  Later  it  was  found,  however,  that  better  devices 

o 

could  be  produced  by  vacuum  deposition  of  first  a  thin  layer  (15A) 

O 

of  chromium  and  then  a  thicker  layer  (1000A)  of  gold.  Vacuum 

deposited  gold  alone  was  not  satisfactory,  since  the  contacts  then 

deteriorated  upon  heating;  apparently  the  chromium  is  needed  to  make 

the  contact  stick  better  to  the  material.  In  this  manner  very 

reproducible  data  could  be  obtained. 

Fig.  2  shows  D*  versus  thickness  w  for  PLZT  samples  with  Cr  -  Au 

-1/2 

contacts.  Approximately  a  D*  =  const  w  relationship  was  found 
in  agreement  with  the  elementary  theory.  The  measurements  were 
performed  at  100  Hz. 

Table  1  gives  D*  measurements  for  PLZT  samples  different  types 
of  contact.  As  the  table  indicates  the  samples  with  Cr  -  Au  contacts 
were  the  best;  the  other  contacts  gave  less  reproducible 
and  had  in  general  a  lower  D*. 
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This  report  summarizes  work  done  on  the  theory  of  ferroelectric  detectors,  sources 
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capacitive  bolometers.  In  addition,  results  of  other  work  including  theory  of  thenral 
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tion,  problems  encountered  in  device  fabrication  and  measurements  of  D*  at  100H* 
of  PL2T  samples  are  summarized. 

The  main  results  are  that  the  parameters  of  the  thermodynamical  theory  ol  pyro¬ 
electric  detectors  can  be  expressed  in  terms  of  independent  f undomental  molecular 
parameters  which  yield  a  lossless  material  figure  of  merit  p/v~ t.  in  close  ugivo- 
ment  with  experimental  valves.  This  approximate  model  loased  on  the  moLecul.ii  field 
approximation  predicts  that  p/<~~£  is  almost  independent  of  the  material  parameters 
over  a  wide  range  of  pyroelectric  materials.  Further  improvement  in  detectivity 
must  come  from  lowering  the  material  loss  parameter  tan<5.  Theoretical  work  on  the 
noise  sources  in  ferroelectrics  led  to  the  conclusion  that  the  thermal  noise  in 
caused  by  the  partition  noise  in  a  two- level  system.  A  capacitive  bolometer  o]xirating 
in  the  paraelectric  region  near  the  Curie  temperature  exhibited  a  very  low  noise 
equivalent  power  which  is  comparable  with  that  of  a  pyroelectric  detector.  It  is 
concluded  that  an  understanding  of  the  fundamental  loss  mechanisms  in  ferroelectric  ' 
pyroelectric  detectors  is  needed  before  improvement  in  a  mnterials-1 imited  device 
detectivity  can  be  attained. 
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